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Wake-driven Dynamics of Finite-sized Buoyant Spheres in Turbulence 
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Particles suspended in turbulent flows are affected by the turbulence and at the same time act back 
on the flow. The resulting coupling can give rise to rich variability in their dynamics. Here we report 
experimental results from an investigation of finite-sized buoyant spheres in turbulence. We find 
that even a marginal reduction in the particle’s density from that of the fluid can result in strong 
modification of its dynamics. In contrast to classical spatial filtering arguments and predictions 
of particle models, we find that the particle acceleration variance increases with size. We trace 
this reversed trend back to the growing contribution from wake-induced forces, unaccounted for in 
current particle models in turbulence. Our findings highlight the need for improved multi-physics 
based models that account for particle wake effects for a faithful representation of buoyant-sphere 
dynamics in turbulence. 


Particulate suspensions in turbulent flows are found in 
a wide range of natural and industrial settings — typical 
examples include pollutants dispersed in the atmosphere, 
droplet suspensions in clouds, air bubbles and plankton 
distributions in the oceans, and sprays in engine combus¬ 
tion [l]-[5’]. The behavior of a particle in a flow is intri¬ 
cately linked to several quantities such as the particle’s 
size and shape, its density relative to the carrier fluid, 
and the flow conditions among others [fl]. For modeling 
purposes, the equations governing particle motion are of¬ 
ten simplified to the case of a dilute suspension of small 
rigid spheres in a non-uniform flow [3, @. In this frame¬ 
work, particle motion in turbulence is described by three 
fundamental control parameters: the ratio of particle size 
to dissipative length scale (S = d p /r )), the particle-fluid 
density ratio (r = Pp/Pf)i and the Taylor Reynolds num¬ 
ber (Re^) of the carrier flow. 


for flow non-uniformity and spatial filtering at the parti¬ 
cle’s scale through corrections such as the Faxen correc¬ 
tions (FC) to the force terms [2]. 
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In many practical situations, particles have a finite- 
size compared to the dissipative length scales of the flow. 
Experimental studies have addressed the effects of finite 
size mainly by using neutrally buoyant finite-size parti¬ 
cles in homogeneous and isotropic turbulence 0, HEH ■ 
These studies highlighted certain effects of finite size on 
the particle’s statistical properties, namely a decrease in 
acceleration variance, an increase in correlation times and 
a decrease in intermittency in the acceleration PDFs (on 
increasing the particle’s size). All three effects could be 
interpreted through classical inertial range scaling argu¬ 
ments, which propose spatial filtering due to increasing 
particle size as the underlying mechanism. From a model¬ 
ing perspective, these effects were captured by accounting 
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FIG. 1: Parameter space of density ratio (r) vs. size ratio (E) 
for particles in turbulence. Datapoints from literature: A - 
Voth et al. Q; > - Volk et al. 13|; * - Gibert et al. [3; 
o - Brown et al. ljl; o - Qureshi et al. [10]; + - Volk et al. |l5l |; 
V - Mercado et al. [l6j ]: • - Prakash et al. 0; and present ex¬ 
periments: □ - Marginally buoyant, U - Moderately buoyant 
and □ - Very buoyant particles. 

The FC model worked reasonably well for finite-sized 
neutrally buoyant particles 0 El, El, 51 0 in turbu¬ 
lence. This also led to its extension to predict the be¬ 
havior of other classes of particles, namely heavy and 
buoyant particles [2, IH, 53, [Hi ■ Building on these predic¬ 
tions, some generic models have been proposed to predict 
the rms of the acceleration of arbitrary-density finite-size 
particles [20]. These extensions, which practically en¬ 
compass many of the naturally and industrially relevant 
particle-laden turbulent flows (where r ^ 1, H > 1), are 
urgently pending experimental validation. Since fully re¬ 
solved numerical simulations (e.g. Physalis j2l| or Front 
tracking [22j]) are too expensive for high Re^, one is in 
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need of experiments in this regime. However, conduct¬ 
ing experiments with non-neutrally buoyant particles has 
been a challenging task. In zero mean-flow turbulence 
setups pi [23l425| these particles would drift vertically 
past the small fixed measurement volumes, making it ex¬ 
tremely difficult to obtain long particle trajectories in the 
Lagrangian frame. 

In this Letter , we present a novel experimental strat¬ 
egy, wherein a mean-flow may counteract the drift of 
the particles. Our investigation covers the regime of 
finite-sized buoyant spheres (4 < S < 50) in turbu¬ 
lence (see Fig. [[]). The size-ratio (S) is defined as the 
ratio between particle diameter d p and the Kolmogorov 
length scale 77 in the flow. We study three density 
ratios (r ~ 0.92, 0.52 & 0.02), where T is the ratio 
of particle density p p to fluid density pf. According 
to the F value, we call the particles either marginally 
buoyant (r ~ 0.92) or moderately buoyant (r ps 0.52) 
or very buoyant (r 0.02). The Galileo number, 

Ga = y/<?dp(l — T)/V, provides a good estimate of the 

buoyancy force in comparison to viscous force. We cover 
two orders of magnitude variation in Ga (sa 30 — 3000), 
and characterize the Lagrangian dynamics of buoyant 
spheres in turbulence. 

The experiments were performed in the Twente Water 
Tunnel (TWT) facility, in which an active grid generated 
nearly homogeneous and isotropic turbulence in the mea¬ 
surement section EHI. The water tunnel was config¬ 
ured to have downward flow in the measurement section, 
and the Taylor Reynolds number of the flow, Re a, was 
varied from 180 to 300. A small number of rigid buoy¬ 
ant spheres ( 0.8 mm < d p < 10 mm) were dispersed in 
each experiment; the volume fraction <j> was kept suffi¬ 
ciently low such that the interaction between the spheres 
was negligible (</> ~ O(10 -5 )). The closed circuit enabled 
the suspended particles to reappear regularly, and hence, 
sufficient statistics could be obtained (see Fig. 1 in the 
supplementary material (13). 

The particles were imaged using two high speed cam¬ 
eras placed at a 90 degree angle between them. The 
spheres appear as dark circles in the back-lit images and 
their diameters corresponded to around 10 pixels. The 
measurement window size was adjusted to ensure this res¬ 
olution, and the circle-centers were accurately detected 
using the Circular Hough Transform method. A parti¬ 
cle tracking code was used to obtain trajectories of the 
spheres. These were further subjected to smoothing us¬ 
ing a spline based technique [35j , which yielded robust re¬ 
sults across the different experimental conditions. Three 
dimensional trajectories were subsequently obtained by 
matching the particle tracks from the individual cameras 
using a cross-correlation based method. 

We first address the question of how a marginal reduc¬ 
tion in density ratio (T « 0.92) affects the dynamics of 
finite-sized particles in turbulence. We focus on one of 
the horizontal components, as in this direction, there is 
no gravity-induced directional preference. In Fig. Ufa), 
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FIG. 2: Acceleration statistics for buoyant particles in tur¬ 
bulence at Re> fs 180. (a) Acceleration PDF for marginally 
buoyant spheres. Inset shows (a 2 ) 1 ^ 2 (mm/s 2 ) vs. E = dp/ 77 . 
The dash-dotted line (in blue) with slope ss -1/3 shows the 
prediction based on classical inertial scaling arguments 0 - 
(b) Normalized acceleration variance from experiments (EXP) 
compared to the results from Faxen-corrected (FC) simula¬ 
tions at Re> = 180 [2]. The horizontal dashed and dash- 
dotted lines (in black) are lines marking the tracer particle 
acceleration ((a 2 ) / (a 2 ) = 1) and the upper bound of FC 
simulations ((a 2 ) / (a 2 ) = 9) respectively. Here, a/ is the 
fluid tracer acceleration. 


we show the normalized acceleration probability density 
function (PDF) for marginally buoyant (MB) spheres at 
different S. The ReA is maintained constant (« 180), and 
the particle size ranges from a few Kolmogorov lengths 
to a fraction of the integral scale [H|. At the smallest 
size ratio (S « 4) the PDF has wide tails, while it is 
nearly Gaussian at 3 ~ 50. This trend of narrowing 
PDF tails with increasing 5 was predicted for buoyant 
particles using the FC model |2!, [Tt} . Similar behaviour 
was reported experimentally and also predicted numeri¬ 
cally (FC model) for neutrally buoyant spheres in turbu¬ 
lence @,m. 

While it seems that accounting for flow non-uniformity 
and spatial filtering are sufficient to capture the PDF 
trends, the absolute accelerations reveal something con¬ 
trary to expectation. The inset to Fig. [U(a) shows that 

the rms of particle acceleration ((a 2 ) ) increases with 

S, suggesting that larger particles experience greater 
acceleration fluctuations. This observation is in con¬ 
trast with both classical inertial range scaling predic¬ 
tions ((a 2 ) 1 / 2 ~ S -1 / 3 ) and the predictions of FC 
model (Fig. (Hb)), according to which larger particles 
should experience milder acceleration fluctuations due 
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to spatial filtering effects [9j. While this —1/3 scaling 
has been experimentally validated for neutrally buoy¬ 
ant (NB) spheres in turbulence 0,1 El, surprisingly, 
only a marginal reduction in particle density reverses the 
trend. Furthermore, the accelerations exceed even the 
upper bound of the Faxen model predictions for buoy¬ 
ant particles i.e. (a 2 ) / > 9 ( see marginally buoy¬ 

ant (r « 0.92) case in Fig. [Hb)). Here, a/ is the mea¬ 
sured fluid tracer acceleration, which was calculated from 
particle tracking experiments (conducted separately) us¬ 
ing ss 20 /im diameter fluorescent tracer particles [411] , 

In order to understand the surprisingly large deviations 
from current particle model predictions, we look into the 
temporal response of the marginally buoyant particles. 
The Lagrangian autocorrelation function C a (r) plotted 
for different particle size-ratios are shown in Fig. [31(a) 
& (b). At lower size-ratios (see Fig. [3](a)), the particle 
accelerations decorrelate according to a turbulent spa¬ 
tial filtering based time scale, Td = (dp/e) 1 / 3 . This 

corresponds to the MB particle with the smallest size- 
ratio (5 « 4) in the inset to Fig. [3](a). At larger size- 
ratios (see Fig. ECb), where S ss 16, 25 and 50), the 
accelerations decorrelate according to a vortex shedding 
time-scale, t v = St ^u , instead of r^. Here, Sth is the 
Strouhal number and U r is the measured mean drift ve¬ 
locity of the particle in the turbulent flow. This decorre¬ 
lation behavior is robust across the different Re^, S, and 
r cases (see table|T]), with the larger particles additionally 
displaying strong periodicity. We therefore see a gradual 
transition from the steady drag regime (Ga ~ 30) to a 
regime with growing vortex shedding-induced effects (Ga 
from 225 to 900). 

For a particle in turbulence, there are two effects 
which can lead to a change in (a 2 ) : (i) the con¬ 

tribution from turbulence - this decreases with increase 
in particle size due to spatial filtering effects and (ii) 
the contribution from vortex shedding - this however in¬ 
creases with increasing particle size. Therefore, the total 
( a 2 ) is the combined effect of turbulence-induced and 
vortex-shedding-induced forcing. As 5 increases from 
4 to 50 (see inset to Fig. [21(a)), the unsteady forcing 
due to vortex shedding starts to outweigh the turbulent 
forcing at the particle’s scale. This results in an over- 

I/O 

all increase in (a 2 ) “ with 5 (for a detailed explana¬ 
tion foiMhe observed trend, see the supplementary ma¬ 
terial US). Thus, the agreeing trend of decreasing in- 
termittency in normalised PDFs (between Fig. [2](a) and 
FC model predictions i) is rather a coincidence, since 
the absolute particle accelerations are orders of magni¬ 
tude different. The present findings reveal the role of 
wake-driven forces (vortex shedding) in buoyant particle 
dynamics in turbulence. 

We now study the effect of the density ratio T on 
buoyant sphere acceleration statistics. In Fig. [2](b), 
we show the results of varying particle density at con¬ 
stant H ~ 50. Not surprisingly, the FC model underpre¬ 
dicts the acceleration variance for all three density ratios. 



t/t v 

FIG. 3: (a) Lagrangian acceleration autocorrelation func¬ 
tion C a {r) for marginally buoyant (MB, T S3 0.92) and 
neutrally buoyant (NB, T « 1) spheres with low Ga (< 30), 
along with fully resolved DNS predictions from Ref. [Till 
for NB particles at Rex = 32. (b) C a (r) for MB parti¬ 

cles (T R3 0.92) with moderately large Ga (225 < Ga < 900) 
at Rca ~ 180. 


TABLE I: A summary of the explored parameter space. 
Ga — Galileo number (approximate), T = p P /pf — density ra¬ 
tio, S = d p / 7 1 — size ratio, td — time corresponding to the first 
minima of C a (r), and t„ — vortex shedding timescale based 
on Sth ~ 0 .2. Re\ lies in the range 180—300. Supplemen¬ 
tary movies A, B, C and D show 3D trajectories of marginally 
buoyant, moderately buoyant and very buoyant particles [27]. 


Particle type 

Ga 

r 

H 

1/td (s x ) 

tv /td 

Very buoyant 

3000 

0.02 

30-60 

8.3-8.4 

1.04T0.01 

Moderately 

buoyant 

2000 

0.52 

30-60 

5.16-5.26 

0.95±0.01 


900 

0.92 

30-60 

2.45-2.69 

1.02T0.05 

Marginally 

450 

0.92 

18-36 

2.5-2.7 

1.04T0.04 

buoyant 

225 

0.92 

12-24 

2.4-2.5 

0.99T0.03 


30 

0.92 

4-8 

4.0-8.0 

3.00T1.00 


The deviation is greatest for the ‘very buoyant’ particle 
(r « 0.02; Ga ~ 3000), whose experimental result is al¬ 
most three orders of magnitude higher, suggesting that 
the particle’s accelerations do not originate from the tur¬ 
bulent forcing. In Fig. |U(a) & (b), we present the hori¬ 
zontal and vertical components of the velocity and accel¬ 
eration PDFs for this very buoyant particle at three Rey. 
The horizontal velocity PDF shows a symmetric flat-head 
distribution, and the horizontal acceleration PDF has a 
bumped-head distribution. Both the velocity and accel¬ 
eration statistics show no observable effect of changes in 
Re\ or B. The vertical components have similar char- 
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FIG. 4: PDFs of (a) horizontal and (b) vertical components of velocity and acceleration for a very buoyant (r « 0.02) finite¬ 
sized sphere at Ga « 3000. Blue, green, and red colours correspond to Rex ~ 180, 250, and 300, respectively. The variable 
X corresponds to Vh, ah, v y and a y in their respective PDFs. The PDFs have been shifted vertically for clarity in viewing. A 
Gaussian profile (thin black line) is overlaid for comparison, (c)-(f) show time traces of velocity and acceleration components 
of a representative trajectory at Rex ~ 300. (g)-(j) show the trajectory corresponding to the time traces in (c)-(f), projected 
on a vertical plane and color-coded with instantaneous quantities. 


acteristics, but the distributions show positive skewness, 
particularly for acceleration. The Lagrangian autocor¬ 
relations of velocity and acceleration are oscillatory (see 
supplementary material S3), with a period comparable 
to t v (see tableUJ). These observations indicate that buoy¬ 
ant particle dynamics is strongly dominated by its wake- 
induced forces. However, the role of turbulence is still 
evident in the PDF-tails. 

In Fig. QJc)-(f), we plot the time-series of velocity and 
accelerations of a representative particle trajectory of the 
very buoyant sphere from our experiments. A careful ex¬ 
amination reveals that the non-regular PDF shapes have 
their origin in the time history of the cyclic motions ex¬ 
perienced by the particle. In particular, the flat-headed, 
velocity, the bumped-head acceleration, and the positively 
skewed vertical PDFs can be traced back to the specific 
details of the periodic motions. In Fig. SDg)-(j), we plot 
the trajectory coresponding to this time trace projected 
on a vertical plane. The particle traces regular Lissajous 
orbits, only weakly disturbed by the surrounding turbu¬ 
lence (see supplementary movie D ^3)- While such peri¬ 
odic motions are classically observed for tethered spheres 
in laminar flows SHI, it is intriguing that these dom¬ 
inate even under highly turbulent conditions, with the 
turbulence playing only a secondary role. 


In summary, the dynamics of finite-sized buoyant par¬ 
ticles in turbulence can be strongly two-way coupled. 
We show that, despite accounting for finite-size effects, 
current particle models break down in situations where 
the particle density is not perfectly matched with the 
fluid. Wake effects are evident in the experiments, and 
particle models which ignore these influences are un¬ 
able to capture the statistics accurately. Traditional ap¬ 
proaches have studied particle dynamics by considering 
wake-induced [39j] and turbulence-induced effects Q as 
separate fields. We emphasize the need to merge these 
two fields to develop realistic multi-physics based La¬ 
grangian particle models. Our results provide the first 
key experimental reference for such emerging models (4Cj j 
in the particles in turbulence community. 
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